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Description of nuclear octupole and quadrupole deformation close to the axial 
symmetry: Critical-point behavior of 224 Ra and 224 Th 

P.G. BizzetE] and A.M. Bizzeti-Sona 
Dipartimento di Fisica, Universitd di Firenze 
I.N.F.N., Sezione di Firenze 
Via G. Sansone 1, 50019 Sesto Fiorentino (Firenze), Italy 
(Dated: February 6, 2008) 

The model, introduced in a previous paper, for the description of the octupole and quadrupole 
degrees of freedom in conditions close to the axial symmetry, is applied to situations of shape phase 
transitions where the quadrupole amplitude can reach zero. The transitional nuclei 224 ' 226 Ra and 
224 Th are discussed in the frame of this model. Their level schemes can be reasonably accounted 
for assuming a square-well potential in two dimensions. Electromagnetic transition amplitudes are 
also evaluated and compared with existing experimental data. 



PACS numbers: 21.60.Ev 



I. INTRODUCTION 



The phase transition between spherical and axially de- 
formed quadrupole shape of nuclei has been the object 
of several theoretical and experimental works in recent 
years. In particular, the properties of nuclei close to 
the critical point, predicted by Iachello's model of X(5) 
symmetry 111, have been actually observed in several 
cases 0, H, 13, [E @, > while some other nuclides show- 
ing the ratio E(4+)/E(2+) w 2.91 expected for the X(5) 
symmetry are presently under investigation. Moreover, 
in the Ra - Th region, it has been observed that the 
isotopes 224 Ra and 224 Th have a positive-parity ground- 
state band with a sequence of level energies very close to 
the X(5) predictions [1, Q. Here, however, the presence 
of a very low lying negative parity band, soon merging 
with the positive-parity one for J > 5, proves that the 
octupole mode of deformation plays an important role 
and should not be ignored in discussing the behavior of 
the phase transition. 

In a previous paper (To| (henceforth referred to as 
I) a simple model has been introduced to describe the 
phase transitions in nuclear shape involving the octupole 
mode 1 . To this purpose, a new parametrization of the 
collective coordinates describing the nuclear quadrupole 
and octupole deformation has been introduced and dis- 
cussed. The nuclear shape is represented in the intrinsic 
frame defined by the principal axes of the overall ten- 
sor of inertia, in situations close to (but not necessarily 
coincident with) the axial-symmetry limit. In the same 
paper, a specific model is developed to describe the crit- 
ical point of the phase transition in the octupole mode, 
between harmonic oscillations and permanent asymmet- 



*Electronic address: bizzeti@fi.infn.it 

1 We have now the occasion to correct a few misprints which es- 
caped proofs revision in paper I: 

Eq. (23c) should read qs = (L3 — p v — 2p x — 3p$) / (4«q) . 
In the Table VII, the 4th element of the 5th line should be 
sin 82 sin 63 / 3l • We apologize for these errors. 



ric deformation, in nuclei which already possess a stable 
quadrupole deformation. The Thorium isotopic chain 
was investigated and the experimental data concerning 
226,228 Th were compared with the model predictions [Icj |. 
The former appears to be close to the critical point, while 
the latter can be interpreted as an example of harmonic 
oscillations in the axial octupole mode. 

In the present paper, we extend the investigation to the 
cases where the quadrupole deformation is not steady 
but performs oscillations under the effect of a proper 
potential, and in particular for situations close to the 
quadrupole critical point described by the X(5) symme- 
try, in the Radium and Thorium isotopic chain. 

As we shall see, the properties of the already men- 
tioned nuclei 224 Th and 224 Ra result to be reasonably 
described by our model with a "critical" (flat) potential 
well, extending both in the /?2 and (3s directions. More- 
over, we observe that, as far as the level scheme is con- 
cerned, also the next isotope 226 Ra can be accounted for 
with a proper critical-point potential, in spite of the fact 
that the positive-parity part of the ground-state band 
does not follow the X(5) predictions. As in the case of 
Thorium, heavier isotopes of Radium have a permanent 
quadrupole deformation and octupole excitations of vi- 
brational character, while the lighter ones are either non 
collective or vibrational in the quadrupole mode. 

Some results of this work, at different phases of ad- 
vancement, have been reported at several Conferences or 

schools dGH, [13, El- 

For convenience of the reader, we report in the next 
Section fll] some evidence of the phase transitions in the 
Radium and Thorium isotopic chain, while the definition 
of variables introduced in I, and a few results relevant 
to the present work, are briefly summarized in the Sec- 
tion |TTT21 In the following subsections, the model intro- 
duced in I is specialized to a form suitable for a critical 
potential in two dimensions. Finally, in the Section IIVI 
the model results are reported and compared with the 
existing experimental evidence for 224,226 Ra and 224 Th. 

Previous models of quadrupole-octupole deformation 
are quoted in I. Since then, new relevant papers have 
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appeared. A new "Analytic Quadrupole-Octupole axi- 
ally symmetric model" (AQOA) has been proposed by 
D.Bonatsos el al. 14] to discuss the evolution of the 
quadrupole and octupole collectivity in Ra and Th iso- 
topes. A parameter free model starting from a simi- 
lar approach has been developed by Lenis and Bonat- 
sos [15| and compared with the experimental results for 
226 Ra and 226 Th. Moreover, a variant of the AQOA 
model, introducing a renormalization of the nuclear mo- 
ment of inertia [161 ]. has been used to describe the low- 
est quadrupole and octupole bands of the N = 90 iso- 
tones 150 Nd, 152 Sm, 154 Gd and 156 Dy. A discussion of 
the octupole bands of 150 Nd and 152 Sm is also contained 
in the Ref.s 0, [12]. Finally^ an extension of the Ex- 



tended Coherent State model 17 1 has been developed by 



IE Hi 



20l | . to include in 
1 + and 1~ bands, 



A. A. Raduta and coworkers 
the model space also the lowest K 71 
and the model predictions have been compared with the 
experimental data for several nuclei of the regions of the 
actinides and of the rare-earths. 



II. PHASE TRANSITIONS IN THE Ra 
REGION 
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FIG. 1: (Color online) Indicators of the quadrupole collectiv- 
ity (left) and of the octupole collectivity (right), as a function 
of the neutron number N in the isotopic chain of Ra (circles) 
and Th (triangles): (a) - Excitation energy of the first 2 + 
level; (b) - Energy ratio E(A + ) / E(2 + ); (c) - Excitation en- 
ergy of the first level of the K n = 0~ band, Jq = 1~ (open 
symbols) and of the lowest known level of other negative- 
parity bands, Jo = 2~ or 1^ (full symbols); (d) - Energy ratio 
E(l~)/E(2 + ). The horizontal line in the part (6) shows the 
value (2.91) expected for the X(5) symmetry (from ref. [l(il]h 



We summarize here the existing evidence for the evolu- 
tion of nuclear shapes for Radium and Thorium isotopes 
in the transitional region N = 130 — 140. 

The Fig. [IJ taken from I, shows the behavior of some 
indicators of quadrupole and octupole collectivity, as a 
function of the neutron number N, in the isotopic chain 
of Ra and Th. It has been noted in I that 226 Th appears 
to be close to the critical point in the octupole defor- 
mation, while it possesses a stable quadrupole deforma- 
tion @2- At larger values of N, Th isotopes maintain a 
stable quadrupole deformation, while the octupole mode 
evolves towards the vibrational behavior, as indicated by 
the large excitation energies of all negative-parity lev- 
els. At N — 130 or less, the quadrupole mode has a 
vibrational (or non collective) character. It turns out, 
therefore, that the octupole phase transition proceeds in 
the direction opposite to the one of quadrupole. We also 
observe that the phase transition only involves the ax- 
ial octupole mode. In fact, the energy of the J n = 1~ 
band head of the K" = 0~ octupole band shows a sharp 
decrease, both in its absolute value and in the ratio to 
the E(2 + ), when the neutron number decreases below 
N = 142. Other octupole bands (with K > 0) do 
not show a similar trend (fig. [Uc)), and one can con- 
clude that non-axial octupole excitations maintain a vi- 
brational character. A similar trend is apparent also for 
Ra isotopes. 

In order to describe Th and Ra isotopes with A < 226, 
we must allow also f3 2 to vary and perform (non neces- 
sarily harmonic) oscillations. If we consider the value 
E(4 + )/E(2 + ) — 2.91 as a signature of the critical point 
with respect to the quadrupole deformation, this would 
correspond approximately to 224 Ra and 224 Th. 



III. THE MODEL FOR QUADRUPOLE 
OCTUPOLE VIBRATIONS 



A. Summary of the variable definitions 

The dynamical variables (A = 2,3; fi = —A. ..A), 
describing the quadrupole and octupole deformation in 
the intrinsic reference frame, are parametrized as 



(2) 



aS 2) = - 



/? 2 COS 72 W #2 



,(2) 



,(3) 



,( 3 ) 



1/2 p 2 sin 72 - i 



v (sin if + i cos ip) 



it sin x 



f3 3 cos 73 w f3 3 

v (sin if + i sin ip) 



,(3) - 



\/l/2 sin 73 + i 



P2 



« sin x 



,(3) 



w sin?? 



+ i w 



cos i? 



cos 73 + (VT5/2) sin 73 
COS73 - (v / 15/2) sin7 3 



w (sin i? + i cos -d) 



(1) 



With this choice, valid in situations close to the axial 
symmetry, the tensor of inertia turns out to be diagonal 
up to the first order in the small quantities describing the 
non-axial deformations. 
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In the Eq.s[T]the variables 72 and 73 are still employed, 
in order to keep some transparency with respect to the 
standard expressions used to describe the quadrupole [2l| 
or the octupole deformation alone (2^]. However, it is 
more convenient to substitute them with expressions in- 
volving the variables u, \ and a new variable uq\ neglect- 
ing second-order and higher-order terms, 



72 



73 = 



10ft 



5/31 



u cos x 



5/?l 



u 



(2) 



V2/3 2 



u cos x + 



V5/(^,/9 3 ) , 



■//() 



It is possible to show that a definite value of the angular- 
momentum component X along the intrinsic axis 3 and a 
definite parity can be associated to the degrees of freedom 
corresponding to the variables v, x (or u, <p or w, $ or uq): 
= 1~ (or 2~ or 3~ or 2 + , respectively). This result 
is independent of the form of the function /(ft, ft) of 
Eq. [3] (which, actually, was left undetermined in I). 



TABLE I: The matrix of inertia Q after the introduction of 
the variables uo, v, u, w, <p, \, and $ (see text). Here, J\ — 
J2 = 3($ + 2/3f),and J 3 =4f 2 (f3 2 ,f3 3 ) ul + 2v 2 + 8u 2 + I8w 2 . 
Only the leading terms are shown. Neglected terms are small 
of the first order (or smaller) in the sub-matrix involving only 
P2, P3, uo, v, u, w, qi and qi\ of the third order (or smaller) 
in the sub-matrix involving only <p, %, $ and q 3 ; of the second 
order (or smaller) in the rest of the matrix. 
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B. The kinetic energy operator 

The classical expression of the kinetic energy has the 
form 



(3) 



where £ = tjh,fa,^v,x,%<P,w, l &,qi,q.2,q3) and 3ii I2, 
q3 are the components of the angular velocity along the 
three axes of the intrinsic reference frame. As in I, we 
adopt here the convention of including the inertial coef- 
ficient B\ in our amplitudes aj^\ that therefore would 



correspond to \fB\a]^ in the usual notations of Bohr. 
The matrix elements of Q, approximated to the most rel- 
evant order, are shown in the Table HI The determinant 
of this matrix turns out to be 



G cx (p 2 2 +2[3j) 2 f\p 2 ,(3 3 ) ul v 2 u 2 w 2 
= Goifafc) ul v 2 u 2 w 2 . 



(4) 



The Pauli recipe for the quantization of the classical ki- 
netic energy gives the Schrodinger equation 



d_ 



9* 



n 2 



[E-V(0}* = (5) 



where g 2 = G = Det Q and ^ stays for the ensemble of 
the variables £ re . 

To our present purpose, this general treatment must 
be specialized, {e.g. with a proper choice of the arbi- 
trary function f{P%,^) in the Eq.s[3]), keeping in mind 
a necessary condition: the Schrodinger equation for the 
quadrupole amplitude, when the octupole amplitude is 
constrained to small values by a proper restoring poten- 
tial, must converge to that of Bohr and therefore, at the 



critical-point, to that of the X(5) model 2 . As we shall 
see in the Section IIV A[ this result is obtained with the 
choice 



'(/3 2 2 



-Pirn 



2/31) 



5/31 



from which one obtains 



G (/3 2 ,/3 3 ) 



(/3 2 2 +/3 3 2 ) 2 (/3 2 2 + 2/?ir 
(/3 2 2 + 5/3 3 2 ) 2 



(6) 



(7) 



C. The critical potential in two dimensions 

Possible landscapes of axial quadrupole-octupole de- 
formation in the Thorium region are exemplified in Fig.[2j 
where the potential energy is depicted as a function of 
the deformation parameters /3 2 and /33. Reported val- 
ues have been obtained by Nazarewicz et al. [24j with a 
Wood-Saxon-Bogolyubov cranking calculation. We no- 
tice that fig WLfi) shows a potential minimum which is 
localized around a fixed value in the /3 2 direction, while 
a flat minimum extends over a sizable interval in the di- 
rection /?3. This is just the "critical" potential for the 
shape transition between octupole oscillation and per- 
manent octupole deformation (combined with a fixed 



2 This choice is different from the one adopted in I to describe 
the critical point in the octupole degree of freedom with a con- 
stant quadrupole deformation: in such a case, in fact, the proper 
limit for small octupole amplitudes does not correspond to the 
X(5) but to the Frankfurt model [23| . valid for small-amplitude 
octupole vibrations of a well deformed nucleus. 
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FIG. 2: (color on line) Potential-energy surfaces in the fli - 
/?3 plane for several Th isotopes, as given by Nazarewicz et 
al. [3. 



quadrupole deformation), corresponding to the ng.s^G?) 
andd;/), respectively. The fig.s Wa)\^b) and^c), in- 
stead, show a different kind of shape transition, proceed- 
ing directly from a fixed, reflection-asymmetric deforma- 
tion (fig(Hc)) to quadrupole-octupole vibrations around 
a spherical shape (figHJo)). 



The potential corresponding to the critical point is not 
shown. It should be somewhere between fig. [DJa) and 
fig. W[b). One can try to approximate the critical poten- 
tial, as usual, with a square well, but now the flat bottom 
of the well should extend over a finite distance in /3 2 and 
03, and be symmetric in /3 3 around /3 3 = 0. The shape of 
the borders is obviously relevant to the result. One could 
imagine shapes like those shown in fig. [3] with dashed or 
dotted lines, but their description would involve at least 
two or three free parameters, and the comparison with 
experimental data could be not very significant. We have 
found, however, that good results are obtained also with 
a simple rectangular shape (solid line in fig. [3]), imply- 
ing only one free parameter, b = fl^/fl^ (apart from a 
common factor of scale). 




FIG. 3: (color on line) Possible shapes for a potential well 
simulating the critical-point potential. The potential-energy 
surface of fig. [2j6) is also shown for comparison. 



IV. RESULTS AND COMPARISON WITH 
EXPERIMENTAL DATA 

A. The Energy eigenvalues 

Now, as a first step, we can evaluate, as a function of 
b, the level energies in the ground-state band and deduce 
the best value of the parameter from a comparison with 
experimental results (fig. [3]). To proceed, we must do 
some assumptions on the behavior of axial and non axial 
modes of deformation. We will assume that 

• Our choice of variables corresponds to independent 
degrees of freedom. 

• Non-axial vibrations are confined to their lowest 
stationary state. 

• An approximation similar to that of the X(5) model 
is valid for the differential equations of all non-axial 
amplitudes: i.e., the differential equation in /?2, /3 3 
can be approximately decoupled from those con- 
cerning the other degrees of freedom. 

Therefore, the complete wavefunction '5 of Eq. [8] can be 
factorized in three parts, as in Eq. 30 of I: 



* = * (/3 2 ,/3 3 ) *! Y JM (Cl) 



(8) 



where the function 'J'i depends on the deformation vari- 
ables different from 02, 

From the Eq. 2] we know that also the determinant 
G is factorized in the same way. Then, the differential 
equation for /3 2 , @3 takes the form 



-1/2 



d 

d/3 2 



G 



1/2 



d 



d 



dfcj d(3 3 



G 



1/2 



d(3 3 



(9) 



This equation can be somewhat simplified with the sub- 
stitution 



*o(/?2,/fe) = <r 1/2 $(/3 2 ,/3 3 ) 
where g oc G Q , to obtain 



(10) 



d 



d 2 



J(J + i) 



dPl + dPl +e Vip2 ^ 3) 3($ + 2$) 
+ V s (p 2 ,fo)) $(/? 2 ,/3 3 ) = 



with 



1 
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~d 2 g d 2 g 








2ff 


[dp 2 + d/3l 



(11) 



(12) 



With the choice of /(/3 2 ,/3 3 ) given in the Eq. [5J from 
Eq. [7] one obtains 



9 oc 



(/5| + 5/5|) 



(13) 
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FIG. 4: (Color online) Experimental excitation energies of the positive parity levels (circles) and of the negative parity ones 
(triangles), in units of E(2 + ), for 224,226 Ra and 224 Th, compared with the results of the present model (full line) with the 
following values of the parameter b = //3™ : 0.81 for 224 Ra, 0.68 for 226 Ra, and 0.85 for 224 Th. The predictions of the X(5) 
model (dotted lines) and for a rigid reflection-asymmetric rotor (dashed-dotted) are also shown for comparison. 



and, for \fa\ < fa, g oc ft [l+4(/3 3 //3 2 ) 4 + ...] . Therefore, 
the first and second derivative of g with respect to fa 
tend to zero when \ fa\ <C fa and, at the limit fa — > 0, 
V g = — 2 as in the original Bohr model. With the substi- 
tution = d 1 ^ 2 ^, an d assuming V(fa,fa) = inside 
the potential well and = +oo outside, the differential 
equation to be solved takes the form 



finite differences are used: 



e.g.. 



9 



9 



dp? 



d(3? 



<*>(fa,fa) = (14) 



with V s given in the Eq. [12] and $ = on the contour of 
the potential well. The numerical integration has been 
performed with the finite difference method. Namely, the 
space is discretized on a rectangular lattice and values of 
<I> at the lattice centers are taken as independent vari- 
ables. In the place of second derivatives, the ratios of 




FIG. 5: (Color online) Calculated energies of excited levels 
of the ground and first excited band, in units of E{%[ ), as a 
function of the ratio b = ffi/ffi!. 



<9 2 $ 



$(x + A x , y) — 2 y) + $(s - A a , y) 



x.y 



As$(fa,fa) = (-l) J $(fa,-fa), it is enough to consider 
only the region fa > 0. The lattice centers are chosen 
as fa = k 2 A x , fa = (k 3 - l/2)Aj, with k 2 = l...n 2 , 
k 3 = I...713, and A x = ^/(n 2 + l), A y = 2^/(2n 3 + l). 
The integration region is the upper rectangle with < 
fa < P2, < fa < fig. At the upper and lateral borders 
of the rectangle, the value of the cigcnfunction must be 
zero. 

The boundary conditions at fa = are not specified, 
but to evaluate the approximate derivatives with respect 
to fa it is enough to consider the value of <!> at the line 
of centers immediately below zero, where they are either 
equal or opposite to the corresponding ones at fa = A y /2 
according to the even or odd value of J. 

The number of centers internal to the integration re- 
gion - and therefore the number of independent values 
of $ - is now TV = n 2 • ^3, and we obtain a finite dimen- 
sional N x N Hamiltonian matrix. This Hamiltonian has 
been diagonalized with the Implicitly Restarted Arnoldi 
- Lanczos method, using the ARPACK package [25| . 

In the fig.0 calculated values of the excitation energies 
(in units of E(2±) ) are depicted as a function of the 
ratio b = At the limit for (3f -> 0, the curves 

corresponding to even J and tt tend to the X(5) values, 
as expected. With increasing 6, at the beginning these 
curves deviate substantially from the X(5) limit, but they 
come closer to the initial values for b sa 1 . In this region 
it is possible to find a good fit of the ground-state band 
of 224 Ra and of 224 Th, for b = 0.81 and for b = 0.85, 
respectively (Fig. [4]) . 

Moreover, a rather good fit of the ground-state band 
of 226 Ra is obtained with b — 0.68, i.e. close to the 
maximum of the curves for even parity and spin. 
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We can observe that, with our choice of the parameter 
b, the calculated 1~ level is always somewhat lower than 
the experimental one (Fig. |4|) . This fact can be related 
to the inclusion, in the potential well, of a region where 
(3s remains large while (3 2 tends to zero. Actually, the 
wavefunction of the first l~ level extends appreciably in 
this region, at variance with other levels of the ground- 
state band. 

It would be of great interest, of course, to extend the 
comparison to the lowest excited band with K — (the 
s = 2 band in the X(5) model notations). Unfortunately, 
in 224 Th no excited + level is known. The non-yrast 
level schemes of 224 ' 226 Ra will be discussed in the follow- 
ing Section [TVD] 



B. Electromagnetic transition probabilities 

Another important test for the model is provided by 
the E2 transition probabilities. The available experi- 
mental information on _B(E2) values is scarce (only two 
transitions in 224 Ra and 226 Ra, one in 224 Th), but we 
hope our work can stimulate interest for new experi- 
mental investigations. The reduced matrix element of 
the quadrupole transition operator A4(E2) between the 
states \s, K = 0, J) and \s', K — 0, J 7 ) can be evaluated 
as 

( S J||M(E2) I s'J') = C 2 (sJ\p 2 \s'J') (J\\Y 2 \\J f ) , (15) 

with (s, J\(3 2 \s', J') = J ^ s J02^s',J' dT and C 2 constant. 
The volume element dr, in our non-cartesian coordinates, 
is the product of the differentials of the coordinate vari- 
ables multiplied by g = G 1 / 2 , with G the determinant of 
the matrix of inertia Q. In our assumptions, the inte- 
grals over all variables apart from (3 2 ad /?3 are inde- 
pendent from one another and from the integral over 
d/3 2 d/33, and their result is 1 (if the corresponding wave- 
functions are properly normalized). As the electric dipole 
and quadrupole operators do not contain derivatives, we 
can exploit the substitution defined in Eq. [10] to express 
the remaining integral as 



0" rPz 

V.jfrVjj' dr = / dfo d/3 3 $ s j/3 2 $ S 'j' (16) 
o 



This integral has been evaluated numerically for values 
of J < 18, with J' = J - 2 (and also with J' = J - 1). 
The reduced matrix element over the angular coordinates 
has the form 



(J\\Yl\\J') = 



(17) 



(-1) J (4^)- 1 / 2 V(2J+1)(2L + 1)(2J'+1) ( 



J L J' 




The absolute values of the ratios of E2 reduced matrix 
elements, Rj{E2) = (J\\M(E2)\\J-2)/(2 + \\M(E2)\\0 + ), 
for transitions within the positive- and the negative- 
parity part of the ground-state band, are depicted, as 
a function of b = (3f / (3 2 , in the fig. [6](a) . Their limit at 
Pf / 02 ~~ * corresponds, as expected, to the X(5) value. 

In addition to the in-band E2 transition, we have to 
consider the El transitions between levels of opposite 
parity. How to treat El transitions in the frame of the ge- 
ometrical model is a big problem, as all El transition mo- 
ments should vanish for a homogeneous fluid of constant 
charge density. In this sense, El transitions are outside 
the Bohr geometrical model. It is usual to assume a con- 
stant electric polarizability of the nuclear matter (27l. [28| 
to obtain the El operator in the form 



(18) 



Finally, the reduced transition probabilities from J to J' 
are obtained as 

£(E2,sJ-> s'J') = (2J+ l)- 1 (sJ||>l(E2)||s'J') 2 . 



This ansatz should be validated by proper microscopic 
calculations. 

Actually, such a calculation has been performed by 
Tsvenkov et al. [29j ] for a number of Radium, Tho- 
rium and Uranium isotopes, in the frame of the Skirme- 
Hartree-Foch model. The electric dipole moment turns 
out to be almost independent of the angular frequency 
in a given isotope, but can change drastically (even in 
the sign) along the isotopic chain. The small value of the 
electric dipole moments in 224 Ra is correctly predicted 
by these calculations. 

Values of the ratios of the El matrix elements, 
Rj(El) = (J\\M(E1)\\J - 1)/(1-\\M(E1)\\0+), obtained 
with the standard form (Eq. [T8f of the El operator, are 
shown in the fig. \@[b) . They reach a maximum for (3™ / (3 2 
somewhat below 1, i.e. just in a region including the val- 
ues assumed for 224 Ra and 224 Th (0.81 and 0.85, respec- 
tively). The calculated values of (Ji\M{EL)\\J f ) for E2 
and El transitions in the ground state bands of 224 Ra 
(6 = 0.81), 226 Ra (b = 0.68) and 224 Th (b = 0.85) are 
given in the upper part of Table [TXJ Values for the cor- 
responding intra-band transitions are very similar in the 
three cases, while the difference can be larger for the weak 
inter-band transitions, as shown in the lower part of the 
Table. 



C. Comparison with experimental transition 
probabilities 

The Table Hill shows a few values of the ratio of reduced 
matrix elements for transitions of the same multipolar- 
ity, that can be deduced from the available experimen- 
tal information. In the same Table, the corresponding 
values calculated with the present model are also shown 
(columns "Crit."), together with the ones expected for a 
reflection-asymmetric rigid rotor ("Rot."). 

The most direct check of the model predictions would 
come from the ratios of B(E2) values in the ground- 
state band. This is possible only in 224 Ra, and only 
for the decays of the lowest 2 + and 4 + levels. With 
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FIG. 6: (Color online) Absolute values of the ratios of reduced matrix elements of the electromagnetic transition operators, 
R(EL, Ji) = M(EL, Ji -> J f ))/ M(EL, L -> 0), with J f = Ji - L, as a function of the parameter b = pf/ffi- Part (a): E2 
transitions; solid lines: even J+ — > J+, starting with 2+ — > + (from the bottom); dashed lines: odd J~ — » J7^, starting with 
3" — > 1". Part 6: El transitions; solid lines: odd J" — » J+, starting with 1" — > 0+ (from the bottom); dashed lines: even 
J+ — » J^T, starting with 2+ — > 1". The vertical lines correspond to the adopted values of the parameter for 226 Ra, 224 Ra and 
224 Th (6 = 0.68, 0.81 and 0.85, respectively). 



the experimental values reported in the NNDC tabu- 
lation 26], B(E2,2+ -> 0+) = 97 ± 3 W.u. and 
B(E2, 4+ -> 2+) = 138 ± 8 W.u., the experimental value 
of the ratio is 1.42 ± 0.09, to be compared with the value 
1.41 that is obtained from the calculated matrix elements 
of Tables HI and Hn] (for b = 0.81). We remind that, in 



TABLE II: Calculated values of the reduced matrix elements 
of El and E2 transitions in 224,226 Ra and 224 Th, normalized to 
those of the lowest lying transition of the same multipolarity. 



(J\M(E1)\J') (J\M(E2)\J') 



Trans. 


224 Ra 


226 Ra 


224 Th 


Trans. 


224 Ra 


226 Ra 


224 Th 


1" 


^0+ 


100 


100 


100 


2+ 


^0+ 


100 


100 


100 


2+ 


<->l- 


149 


147 


150 


3" 


^1" 


127 


129 


127 


3" 


^2+ 


187 


184 


187 


4+ 


^2+ 


164 


166 


163 


4+ 


^3" 


238 


228 


241 


5" 


^3" 


179 


182 


178 


5" 


^4+ 


273 


262 


276 


6+ 


^4+ 


211 


217 


209 


6+ 


^5" 


333 


311 


338 


7" 


^5" 


223 


230 


221 


7" 


^6+ 


371 


346 


376 


8+ 


^6+ 


252 


264 


248 


8+ 


^7" 


435 


402 


439 


9" 


^7" 


264 


274 


261 


9" 


^8+ 


476 


438 


480 


10+ 


^8+ 


289 


304 


284 


10+ 


^9" 


538 


500 


539 


11" 


^9" 


303 


315 


299 


11" 


^10+ 


582 


538 


582 


12+ 


^10+ 


325 


342 


320 


12+ 


<->ll" 


640 


604 


635 


13" 


<->ll _ 


340 


354 


336 


13" 


^12+ 


685 


644 


680 


14+ 


^12+ 


360 


377 


354 


14+ 


^13" 


738 


709 


730 


15" 


^13" 


376 


390 


371 


15" 


^14+ 


783 


750 


775 


16+ 


^14+ 


393 


411 


388 


16+ 


^15" 


833 


812 


822 


17" 


^15" 


410 


424 


405 


17" 


^16+ 


878 


854 


867 


18+ 


^16+ 


426 


443 


421 


18+ 


<->17 _ 


925 


912 


912 


19" 


<-»17" 


442 


457 


437 


o + 


~ir 


84 


42 


150 


o 2 " 


^2+ 


8 


24 


6 




^0+ 


31 


31 


31 




~17 


38 


37 


38 




^2+ 


49 


50 


47 


1J 


~37 


42 


43 


43 






22 


63 


21 


2+ 

Z 2 


-0+ 


15 


17 


16 


2+ 


~ir 


113 


86 


111 


2+ 

Z 2 


^2+ 


18 


25 


17 


2+ 

Z 2 


«->3j 


143 


103 


142 


2+ 

Z 2 


^4+ 


8 


32 


4 


2+ 

Z 2 




33 


71 


34 


2+ 

Z 2 


^0+ 


91 


82 


92 



the X(5) model [l|, this ratio would be 1.59. For 226 Ra, 
the lifetime of the 4+ state is known, but for the first 
excited state only an approximate value (without error 
estimate) is reported. Also in this case, the deduced ratio 
is consistent with the theoretical estimate (see Table ITTTT) . 
These results is encouraging, but would obviously need to 
be validate by a more extensive check, involving higher- 
lying levels, which, at the moment, is not possible. 

A comparison of the two El transition from the lowest 
level l - to the + and to the 2+ states is possible for 
the three isotopes, as well as for the El branches from 
the 3" in 226 Ra. All these amplitude ratios for transi- 
tions within the ground-state band, shown in the upper 
part of Table Mil are in very good agreement with the 
calculated values. However, they are not significantly dif- 
ferent from those expected for a rigid asymmetric rotor 
(as shown in the last column of Table IIIip and also from 
those reported by Lenis and Bonatsos [15[ on the basis of 
a rather different model. We note that, when the tran- 
sitions to be compared have the same multipolarity, the 
model predictions are parameter free, or - more exactly 
- only involve the model parameter f3™ /flf. 

Instead, when the comparison concerns the ratios of 
the reduced matrix elements for El and E2 transitions 
deexciting the same level, the model predictions include 
a further normalization factor (the ratio of constants C\ 
and C2 of Eq.s ll5ll5|) which needs to be determined from 
the experimental data. This comparison is therefore less 
direct, but it is perhaps more significant, as we shall see 
in the following. 

Results concerning the E1/E2 branches in the ground- 
state band are shown in the Table IIVI and also depicted 
in Fig. [51 Experimental values of El /E2 branching ratios 
in 224 ' 226 Ra and 224 Th include those given in the NNDC 
tabulation (2(| [3l| and later results from Ref. [30] . From 
these branching ratios we have deduced the absolute ra- 
tios - given in the "Exp." columns of Table IIVI - of 
the reduced matrix elements of El and E2 transitions, 
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TABLE III: Experimental and calculated values of the ratios of reduced amplitudes of two El or two E2 transitions. 















( Ja|A*(EL)| J'a)/ 


{ Jb 


II \/t (~FT \ \\ V \ 
\M\JbU)\ J B ) 










Transitions 




224 Ra 




226 tj 

Ka 




224 Th 








Ja 


-Jl 


Jb 


- J B 


Experim. 


Crit. 


Experim. 


Crit. 


Experim. 


Crit. 


Rot. 


El 


K 


-2+ 




-o+ 


1.52 ±0.14 


1.50 


1.36 ± 0.12 


1.47 


1.49 ±0.26 


1.50 


1.42 


El 


3r 


-4+ 


sr 


-2+ 






1.11 ± 0.18 


1.24 






1.15 


E2 


4 + 


-2+ 


2f 




1.60 ±0.05 


1.63 


« 1.76 


1.66 






1.60 


E2 


1 2 " 


-37 


1 2 " 


-»ir 


0.71 ±0.10 


1.10 












El 


1 2 " 


-2+ 


1 2 " 




1.49 ±0.16 


1.57 


1.24 ±0.09 


1.62 








El 


2+ 

Z 2 




2+ 

Z 2 








1.29 ±0.08 


1.20 









TABLE IV: Experimental and calculated values of the ratios of reduced amplitudes of El and E2 transitions from the same 
level (in units of their Weisskopf estimates). The columns of calculated values are normalized to obtain the best fit to the 
experimental values for the transitions within the ground-state band. 



[(M(El)/Mw(El))/(M(E2)/M w (E2))] x 10 3 



Transitions 




224 


Ra 




226 


Ra 




224 


Th 




El 


E2 


Experim. 


Crit. 


Rot. 


Experim. 


Crit. 


Rot. 


Experim. 


Crit. 


Rot. 


37-2+ 


37- 


ir 


0.69 ±0.14 


0.55 


0.57 














57-4+ 


57- 


37 


0.98±0.26 a 


0.63 


0.60 


1.36 ± 0.23 6 


2.13 


2.42 








6+-57 


6+- 


4 + 














7.98 ± 1.17 


7.18 


8.21 


77-6+ 


7r^ 


•57 


0.56 ±0.09 


0.66 


0.65 


2.51 ±0.15 


2.23 


2.56 








8+ -77 


8+- 




< 1.22 


0.68 


0.66 








7.19 ±0.72 


7.86 


8.50 


97-8+ 


97- 


77 


< 1.71 


0.71 


0.67 


2.82 ±0.34 


2.38 


2.67 








10+ -97 


10+ — 


8i 














7.78 ± 0.43 


8.41 


8.67 


117-10+ 


117- 


97 








2.76 ±0.27 


2.53 


2.67 


9.35 ± 0.62 


8.64 


8.73 


12+ -117 


12+- 


10+ 








2.85 ±0.25 


2.62 


2.68 


9.06 ± 0.47 


8.83 


8.78 


137 — 12+ 


137- 


nr 








2.15 ±0.29 


2.70 


2.70 


8.45 ± 0.42 


9.01 


8.82 


14+ -137 


14+ — 


12+ 








2.58 ±0.17 


2.79 


2.71 


9.84 ±0.51 


9.15 


8.86 


157 -14+ 


157 — 


137 








2.53 ±0.17 


2.94 


2.73 


9.69 ±0.65 


9.29 


8.89 


177 — 16+ 


177- 


157 








2.78 ± 0.43 


2.98 


2.74 


10.47 ± 1.34 


9.52 


8.92 


18+ -177 


18+- 


18+ 








3.22 ±0.21 


3.05 


2.74 








X 2 1 n (with 


n = 8 


) 










2.13 


1.34 




1.17 


2.03 


Confidence level (%) 










<5 


18 




31 


<5 



"from NNDC [Slon lv. The 5" — 3~ (142 keV) 7 ray observed in 
the reaction data |30{ appears to be contaminated by a close-lying 
transition from a different reaction, as it results from the intensity 
mismatch in the 5~ — » 3 _ — + 2 + (-^ 1 _ ) cascade. 

b from Ref [H. Not included in the fit. 



each of which is expressed in units of the correspond- 
ing Weisskopf estimate, M W (EL) = (47r)- 1 - /2 [3/(i ± 
3)] (l.2A l / 3 ) L e fm L . In the same table are also shown 
the results of the model calculation at the critical point 
(Crit.), which have been normalized to obtain the best 
fit with the experimental values within the ground-state 
band of each nucleus. Values expected for a rigid asym- 
metric rotor (Rot.), normalized in the same way, are also 
shown. The 226 Ra point at Ji = 5 which, according to the 
authors themselves [3l|, could be considered as a lower 



limit, has not been included in the fits. 

For the ground-state band of 224 Th (Fig. I7I8[) . we find 
a satisfactory agreement between the experimental values 
and the model predictions. In this case we have enough 
data to perform a \ 2 t es t of goodness of fit, and we 
obtain x 2 /N = 1.17 with N — 8 degrees of freedom, 
corresponding to a confidence level of 31%. A fit with 
the rigid-rotor values would give a much larger value 
X 2 /N — 2.03, and a confidence level below 5%. Also for 
the ground-state band of 224 Ra, the few available exper- 
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imental values (or limits) are not far from the results of 
the model, but more experimental data would be nec- 
essary for a significant comparison. Actually, as it was 
soon recognized [HI, [H, [34[ the El transitions in 224 Ra 
are rather weak compared to other nuclei in this region, 
and in particular their strengths are two orders of mag- 
nitude smaller than the corresponding ones in 224 Th. 

Instead, experimental values for 226 Ra deviate signif- 
icantly from the model predictions and approach those 
expected for a rigid rotor. This fact, combined with the 
slight upward deviation of level energies from the calcu- 
lated curve for J > 14, suggests that the critical point 
of the phase transition in the Ra isotopic chain can be 
situated somewhere below A = 226, and probably close 
to A = 224. 



D. The first excited K = band 



224 



Th 



16+ 


2398 


14+ 


1959 




22 


12+ 


1550 




22 


10+ 


1174 




26 


8+ 


834 




33 


6+ 


535 




54 


4+ 


284 




100 


2+ 


98 


0+ 


100 „ 







16^ 


2433 


u 15~ 


2165 




28 




31 

1739 


14+ 


1994 






24 




33 


12+ 


1587 


> n 


1347 




23 




23 


10+ 


1212 


V 9 _ 


998 




23 






8+ 


870 


v T 


700 




30 






6+ 


562 




465 




58 








3^ 


306 


4+ 


298 


/ 1- 


251 










2+ 


100 

98 






0+ 


100 „ 





2213 




32 


j 13" 


1793 




30 


vll 


1409 


» 9 


25 

1064 


; 7~ 


17 

762 


5- 


7 

511 


3^ 


1 

320 



■&A\ I' 206 



Experimental 



Model 



FIG. 7: (Color online) Level scheme of 224 Th, compared with 
the model predictions for b = 0.85. Calculated and experi- 
mental branching ratios are reported for each level. Experi- 
mental data are taken from the NNDC tabulation [26j . Ex- 
perimental branching ratios from the 7~ and 9~ levels are 
not known. Theoretical values of the level energies (in keV) 
are normalized to that of the level; those of the branching 
ratios are deduced from the matrix elements of Table [TT] with 
the experimental values of the transition energies. Calculated 
branches lower than 1% are not shown. 



As anticipated in Section IIV Al no experimental in- 
formation is available for non yrast levels of 224 Th. For 
224 ' 226 Ra isotopes, a few non-yrast levels are known from 
f}~ decay of 224 > 226 Fr, from a decay of 228 < 232 Th or from 
the 226 Ra(t,p) reaction. Unique assignments of the spin 
and parity have been reported only for part of them. 
Some of these levels, which could be considered as mem- 
bers of the excited K = band (the s = 2 band, in the 
X(5) expression) are reported, together with those of the 
yrast band, in Fig.[9l where also the main decay branches 
are indicated. In the same figure, the model predicted 
levels, and their expected 7 branches, are also shown. 

We can immediately observe that non yrast levels pre- 
dicted by the model are always lower than the experimen- 
tal ones (but a comparably large discrepancy is observed 
also in the s = 2 band of X(5) nuclei @, [H, 0])- In the 
lower part of Table IIII1 the calculated amplitude ratios 
for transitions from the excited K = band are com- 
pared with the corresponding experimental ones, if the 
levels j" and 1^~ shown in the Fig. [9] are interpreted as 
belonging to it. Only the ratio of the two El transitions 
from the 1% level of 224 Ra and from the 2^ level of 226 Ra 
are well consistent with the calculated value, while the 



i?(El/E2) 

0.010- 



0.005- 




10 



15 20 

Ji 



FIG. 8: (Color online) Ratios of the absolute value of the 
transition matrix elements (normalized to the Weisskopf Unit) 
for El and E2 transitions in the g.s bands of 224,226 Ra and 
224 Th, from Ji to J t -1 and Jj-2, respectively: Rj(El/E2) = 
{J i \[M(El)\J i -l)/Mw(El)]/[(J i \M(E2)\J i -2)/Mw(E2)]. 
The dotted lines join the calculated values of the ratio (nor- 
malized to obtain the best fit with the ensemble of exper- 
imental values). The dashed lines join the values expected 
for a rigid rotor. The corresponding values deduced from the 
parameter free model of Ref. |Tf| are (apart for a possible 
staggering between even and odd Ji) almost identical to the 
rotational ones for large values of Ji (Ji > 7) and, for de- 
creasing values of Ji, their trend reaches a minimum around 
Ji = 6 and then increases slightly at lower values of Ji. 
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15 1794 



1 916 




Experimental Model Experimental Model 

224 Ra 226 Ra 

FIG. 9: (Color online) Partial level schemes of 224 Ra and 226 Ra, with the experimentally observed 7 transitions, compared with 
the results of model calculations (with b — 0.81 and = 0.68, respectively). Theoretical level energies (in keV) are normalized 
to that of the first excited level. Experimental energies for the lower levels of 224 Ra are taken from the NNDC tabulation [26| ]. 
those of the 10 + , 12 + and higher levels are deduced from the 7-ray energies given by Cocks et al. [U- For 226 Ra, those of 
levels up to 5~ are taken from NNDC or ref. [3lJ, those of higher levels from Cocks et al. [jjjj. Gamma branches lower than 5% 
(or reported as upper limits) are shown as dotted lines, those between 5% and 25% as dashed lines. Calculated branches lower 
than 1% are not shown. For a comparison of experimental E1/E2 branches with the model prediction at the critical point, see 
Table HV] and Fig. |8] 



corresponding ratios for the two El transitions from the 
I2 level of 224 Ra and for the two E2 transitions from the 
I2 level of 224 Ra seem to be significantly different from 
the model predictions (although the latter is subject to 
a large uncertainty, due to the presence of a competing 
Ml component in the 1^ — > 1^ transition). 

As for the E1/E2 ratios for inter-band transitions, it is 
not obvious that the value of the parameter C1/C2 ought 
to be the same as for transitions within the ground-state 
band, but if we assume to be so, the E1/E2 ratios in 
the decay of the 1^~ level of 226 Ra differ by a factor of 
2 from the calculated values: the ratios to the E2 am- 
plitude I2 — > 3^, with the normalization used in the 
the Table El are (0.75 ± 0.7) lO" 3 for the 17 -> Of El 
transition and (1.05 ± 0.07)10" 3 for the 17 27, to 
be compared with the theoretical values 0.29 • 10 -3 and 
0.46 • 10~ 3 , respectively. 

Therefore, if the first two levels of the excited s = 2 
band are tentatively identified with the 0^" and 1^~ levels 
of 224 Ra, their properties are not so well accounted for. 
For this fact, one can hypothesize different explanations. 
First, we remark that the identification of these levels as 
members of the band can be put in discussion. Actu- 
ally, the O2 level could result from other (collective or 
non collective) modes of excitation, as, e.g., pairing vi- 
bration [HlHil, while the 17 could correspond to (or be 
mixed with) the band head of the K v = 1~ band. Other- 



wise, the observed disagreement could indicate that our 
model is unable to correctly predict states outside the 
ground-state band, in particular if they are not far from 
levels of the non axial modes having the same J T . The si- 
multaneous investigation of axial and non axial modes, as 
it has been performed, via the Extended Coherent State 
model, in the ref.s [13, El, O H3] , is outside our present 
possibilities. 

V. CONCLUSIONS 

An extension of Iachcllo's X(5) model to the axial 
quadrupole + octupole deformation has been developed 
with the formalism introduced in our previous paper 
I [lGj. Assuming that both fa and fa can vary within a 
two-dimensional well with rectangular borders, and with 
a proper determination of a free function of the model, 
the results are found to converge to those of X(5) when 
the interval available for fa tends to zero. The formal- 
ism is therefore suitable to describe the critical point of 
phase transitions involving at the same time the axial 
quadrupole and octupole deformation. 

As anticipated in I, the principal aim of this second 
part of our work was the description of the transitional 
nuclei 224 Ra and 224 Th, which were proposed to be close 
to such a critical point. 
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Actually, in spite of the admittedly crude schematiza- 
tion of the bidimcnsional potential, the relative values 
of the excitation energies of levels (of positive and neg- 
ative parity) in the ground-state bands of both 224 Ra 
and 224 Th are satisfactorily reproduced by adjusting the 
only available parameter (the aspect ratio b — [}™ / 02 of 
the potential well). A good agreement is obtained with 
b = 0.81 for 224 Ra and with b = 0.85 for 224 Th. More- 
over, a good agreement is also obtained for the first part 
of the ground-state band of 226 Ra with a lower value of 
the parameter, b — 0.68: only above J = 14 the experi- 
mental points deviate slightly from the calculated values, 
in the direction of the rigid-rotor curve (Fig. [4]). 

The (few) known ratios of transition strengths in the 
ground state band, for electromagnetic transition of equal 
multipolarity (either E2 or El) are in agreement with the 
model predictions. Unfortunately, only in a few cases 
the ratio of the reduced strengths for transitions of equal 
multipolarity can be deduced from the experimental data 
(see Table IIII|) and in these cases the values expected at 
the critical point are not very different from those of the 



rotational model. 

In some more cases, the relative strength of two tran- 
sitions of different multipolarity (El and E2), coming 
from the same level, can be deduced from the measured 
branching ratio. The comparison with the model requires 
in this case one more parameter, which has been deter- 
mined by a best-fit procedure (see Table ITVl and Fig. [5]). 
But, in this case, the expected trend at the critical-point 
is significantly different from that of a rigid rotator. 

The calculated critical-point values of the ratios 
E1/E2 are in a rather good agreement with the exper- 
imental results in the case of 224 Th (Fig. [8]), while for 
226 Ra the trend of empirical values is closer to the one 
expected for a rigid rotor. For 224 Ra, the El transitions 
are very weak and experimental data are too scarce to 
permit a significant comparison with the model predic- 
tions. 

New and more extensive measurements of the transi- 
tion strengths either in 224 Ra or 224 Th would be highly 
desirable, for a more significant test of the model. 
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